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We compute the finite size effects in the ground state energy, equivalently the effective cen- 
tral charge c e g, based on S-matrix theories recently conjectured to describe a cyclic regime of the 
Kosterlitz-Thouless renorrnalization group flows. The effective central charge has periodic prop- 
erties consistent with renorrnalization group predictions. Whereas c e ff for the massive case has a 
singularity in the very deep ultra-violet, we argue that the massless version is non-singular and 
periodic on all length scales. 
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I. INTRODUCTION 
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' Over the last few years a novel renorrnalization group (RG) scenario with limit-cycle behavior has been discovered 
^v^j , in various models in a variety of physical contexts, including nuclear physics [lj, quantum field theoryjp, , quantum 

■ mechanics 4j, superconductivity [3, Bose-Einstein condensation Q, and effective low energy QCD [3. The subject 
of duality cascades in supersymmetric gauge theory § is also suggestive of limit-cycle behavior. The possibility of 
limit cycle behavior in the RG flow was considered as early as 1971 by Wilson 0, however at the time no models 
with this behavior were known. The possibility of chaotic flows has also been recently considered P.IToj. 

As is readily realized, models with cyclic RG flows have novel physical properties in comparison with theories with 
fixed points. For theories with bound states, the general signature of the cyclic RG is an infinite tower of bound 
\ states with energies related by discrete scaling relations 0, 0,0,0) i- e - a self-similarity of the spectrum under a 
^-H ' discrete scale transformation. This property was termed Russian doll scaling in |{| for obvious reasons. On the other 
hand, for continuous quantum field theories without bound states, the cyclic RG is manifested in the periodicity of 

■ the S-matrix as a function of energy , as anticipated already by Wilson . 

The present article is a continuation of In the latter paper, an exact S-matrix for the cyclic regime of the 
Kosterlitz-Thouless flows was proposed for the first time. Using the S-matrices proposed there, here we investigate 
O ' the finite-size effects, or equivalently the finite temperature effects, of the ground state energy E(R) on a circle of cir- 
cumference R (temperature 1/-R). Conventionally one defines the effective central charge c e ff as E(R) = —irc c g(R)/6R. 
1-^ ■ The quantity c e ff governs many physical properties, including the specific heat |lllll2| . It is known that for theories 
with fixed points, the quantity c c s tracks the RG flow, namely, as a function of R it smoothly interpolates between 
the ultra-violet (UV), R = 0, and infra-red (IR), R = 00, values of the Virasoro central charge of the fixed points 
|13|. For this reason c e g is an interesting property to study for a model with an RG limit cycle, and this is the main 
subject of this paper. As we will show, it has a periodic structure as a function of log R. 

There is another c-function c z (L) which also tracks the RG flow between fixed points, this being the content of the 
celebrated c-theorem which is a statement of the irreversibility of the RG flow. Specifically the theorem states 

that c z (L) is a monotonically decreasing function of scale L in the flow toward the IR. This by itself seems to rule out 
cyclic RG flows in unitary theories, since straightforward RG arguments show that c z {L) should be periodic in logL. 
In fact the proof of the theorem constructs c z as a function of the couplings <?, and the flow of c z is induced from the 
flow of g: dc z /dlogL — (3 g d g c z (g), where f3 g is the [3 function; thus if the couplings are periodic, so is c 2 . Whereas 
c c s(R) is related to a one-point function of the trace of the stress-energy tensor in a finite geometry (see below), 
c z (L) is related to a two-point function with fields separated by L in infinite volume, and they are thus different. 
Nevertheless, for models with fixed points both c e g(R) and c z (L) have essentially the same behavior, and we expect 
the same to be true for a theory with a cyclic RG. This leaves us with a paradox, since our theory is a perturbation of 
a unitary conformal field theory by a hermitian operator. This issue will not be dealt with any further in this work. 

Let us now summarize our main results. In section II we outline the main features one expects for integrable, 
relativistic quantum field theories in 2d that have a cyclic RG, using only standard RG arguments in a model- 
independent way. There a clear distinction is made between massive and massless scattering theories. In section III 
we review the definition of the models and RG flows, extending the arguments to higher level k current algebras. We 
point out two possible limiting theories, one massive, the other massless, depending on whether one has a UV or IR 
fixed point in the isotropic limit. In section IV we describe the proposed S-matrices. The massive S-matrix is the 
same as considered in |3|, whereas the massless one is new. Both these S-matrices are periodic in rapidity, which 
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as explained in section II, is a clear signature of the cyclic RG. In section V we use thermodynamic Bethe ansatz 
techniques to study the finite size effects for the massive case. Though periodic structures consistent with the RG 
flows are found for c e ff, in the deep UV it develops a singularity whose nature is explored in some detail. We derive an 
approximate analytic expression for c e g in terms of Riemann's zeta function, eq. (|65|) which is surprisingly simple and 
agrees very well with exact numerical results. The finite size-effects for the massless case are discussed in section VI; 
however at the present time the appropriate thermodynamic Bethe ansatz equations are unknown, and our analysis 
is thus incomplete. However on general grounds we argue that in the massless case c Q g should be an exactly periodic 
function of logi? on all scales and has no singularities, and thus appears to be more consistent with the RG analysis. 

II. PHYSICAL CONSEQUENCES OF AN RG LIMIT CYCLE 

In this section we describe some general properties that a model with an RG limit cycle should exhibit, as a way 
of anticipating our subsequent results. 

Generally speaking, cyclic behavior in the RG flow can in principle exist at all scales, or can be approached 
asymptotically in the ultra-violet or infra-red, the latter being UV or IR limit-cycle behavior. Once the flow is in the 
cyclic regime, the couplings are periodic 

g(l + \)=g(l) (l) 

where I = logL and L is the length scale. Above, the period A is fixed and model-dependent. In our model A is an 
RG invariant function of the two running couplings. As we will review in the next section, our models have a periodic 
RG on all scales rather than a limit cycle. 

Consider first a theory of massive particles. In the deep IR, the particles appear infinitely heavy and decouple, 
leaving an empty theory which cannot support a low-energy limit cycle. Therefore, we expect that a massive theory 
can only support a UV limit cycle. A massless theory on the other hand can be non-trivial in the IR and thus support 
an IR limit cycle or a limit cycle on all scales. We will explore both possibilities subsequently. 

A. RG equations for the S-matrix 

It is well-known that the RG leads to certain scaling relations for the correlation functions (Callan-Symanzik 
equations). Since the S-matrix is obtained from momentum space correlation functions, it thus also obeys RG scaling 
relations. 

Let S denote the 2-particle to 2-particle S-matrix. For an integrable quantum field theory in 2 space-time dimensions, 
S only depends on the kinematic variable E% m = (Pi + P2) 2 where Pi,2 are the energy-momentum vectors of the 
incoming particles. (For an integrable theory, the incoming and outgoing momenta are the same.) Standard RG 
arguments applied to the S-matrix lead to the scaling relation: 

S(e- l E cul , g(l )) = S{E cm , g(l + l ) (2) 

(See e.g. ref. |16||.') For a theory with a limit cycle, when I = A the above equation implies a periodicity in energy: 

S(e- x E cm ,g) = S(E cm ,g) (3) 

Let us now specialize to 2d kinematics. First consider the massive case, where the energy momentum can be 
parameterized in terms of a rapidity f3: 

P = mcosh/3, p = msinh/3 (4) 

where m is the mass of the particles. The center of mass energy is 

E% m = 2m 2 (1 + cosh/?), (3 — (3\ — f3 2 (massive case) (5) 

Suppose this theory has a UV limit cycle. At high energies (3 is large and E cr[L me 13 / 2 . The relation eq. Q then 
implies a periodicity in rapidity: 

S{fi-2\) = S[fi) (6) 
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We turn now to the massless case. Here the massless dispersion relations E = ±p can be parameterized as: 

E = ^e^ R 1 P=-^-e l3rt , for right — movers 

E = —e~ /3L , P=— "o" e > for left — movers (7) 

where now m is an energy scale. The center of mass energy for a right-mover with rapidity (3r scattering with a 
left-mover of rapidity (3l is 

£ cm = me^ 2 , I3 — Pr — (3l, (massless case) (8) 

If Srl(P) is the S-matrix for the scattering of right-movers with left- movers, then eq. (0 again implies a periodicity: 

S RL (f3-2\) = S RL ((3) (9) 

One sees then that the main difference between the massive and massive case is that in the massive case the 
periodicity in rapidity eq. © implies a periodicity in energy E cm only for large E cm , whereas in the massless case 
it leads to a periodicity at all energy scales. Thus a massless theory with the periodicity eq. is consistent with a 
cyclic RG flow on all scales, whereas the massive case at best can correspond to a theory with a limit cycle only in 
the deep UV. 



B. RG equations for finite size effects 



Let us consider the quantum field theory on a cylinder of finite circumference R and length L going to infinity. 
Viewing the finite circumference as the space the hamiltonian lives on, and the length L as the time, the partition 
function behaves as 

Z ~ e- LE ^ (10) 

where E(R) is the ground state energy of the hamiltonian at finite size R. This ground state scaling function contains 
a lot of information about the theory. From E{R) one conventionally defines the effective central charge c e g: 

m = (id 

In theories with fixed points, c c s(R) is equal to the Virasoro central charge c at the fixed point, and smoothly 
interpolates between the ultra-violet and infra-red values of c. For this reason, it is very interesting to study c c ff in a 
theory with no fixed points but rather limit cycles. 

In order to derive an RG scaling equation for c c s we first relate it to a correlation function. Let (Tff) denote the 
one-point function of the trace of the stress-energy tensor on the cylinder. It is a function of the couplings g and R, 
and is related to c c s as follows |l3j : 

RG arguments applied to the one-point function give the scaling relation: 

(TP (e l R, g(l )) = e^ 1 (T£) (R, g(h + 0) (13) 
Taking / to be equal to the period A, one obtains: 

c/ eS (e*R) = c> eB (R), < ff =^g| (14) 



i.e. c' oS (R) is a periodic function of logi? with period A. 
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III. THE MODELS AND THE CYCLIC RG FLOWS 



In this section we describe an action for our model and review the beta-function and RG flows. Consider a conformal 
field theory with su(2) symmetry, described formally by the action S c f t - The conformal symmetry promotes the su(2) 
symmetry to decoupled left-right current algebra symmetry |l7j . with left-moving (right-moving) currents J a (z) 
(J (z)), a = +, — ,3 where z = x + it,~z = x — it are euclidean light-cone space-time variables. The currents are 
normalized to satisfy the following operator product expansion: 

^M 3 (°) ~ A, ^(^(0) ~ ±^(0) J + (z)J-(0) ~ ^ + -J 3 (0) (15) 

■ 2 z _ ■„■ z 

where k is the level, and similarly for the right-moving currents. The simplest realization of the current algebra is in 
terms of a pair of Dirac fermions, J a = ip' ! cr a 'ip, where a a are Pauli matrices, and corresponds to a k = 1 reducible 
representation with Virasoro central charge c = 2. Our model is based on the irreducible level k = 1 current algebra 
which can be represented in terms of the left-moving part of a free boson $ = 4>(z) + <f>(z): 

J 3 = -^=d z <j>, J ± = ~}= exp (±iV2(j>) (16) 
v2 v2 V / 

The action S c ft is then just the action for a free massless boson: 

5 cft = ^ J<?x\{d^Y (17) 

and leads to ($(x)$(0)) = — log \x\ 2 . The above conformal field theory corresponds to Virasoro central charge c = 1. 
(In the condensed matter context, the c = 1 theory can be obtained from the c = 2 theory by spin-charge separation.) 
Our model is defined as an anisotropic left-right current-current perturbation of the conformal field theory: 

S = S clt + J (lg±(J + T + J-T) - 4g || J 3 J 3 ) , (18) 

where g±, g\\ are marginal couplings. The RG beta- functions determine how the couplings gu , g± depend on the length 
scale e l . These flows possess an RG invariant Q satisfying dQ/dl = 0. At one-loop, Q = g 2 — g\. In a particular 
renormalization scheme, an all-orders beta-function has been conjectured |l8j |. and the higher order corrections to Q 
are @: 

2 2 

o = 9 jzA (19) 

Q (l + g^d-glY [U) 
The cyclic regime corresponds to Q < 0. Let us parameterize: 

Q = (20) 

with h > 0. The coupling h is the main parameter of the theory. As shown in ref. 0, it governs the period of the 
RG cycles: 

g(X + l)=g(l), X=^- (21) 

where g = g\\.±- In each cycle gu flows to — oo and jumps to +oo. As described in 0] one can view the coupling 
constant — arctan(4<7|i//i) as living on the universal cover of the circle. The physical quantities, such as the S-matrix 
and the finite size effects are all finite and as we will see reveal the RG cycles in a completely smooth manner. The 
above periodicity is not approached asymptotically, but is present on all scales; this will be an important criterion 
when we analyze the finite size effects. 

The limit h — > corresponds to being on one of either of the two separatrices <?|| = ±g±, which can be taken in the 
half-plane g± > due to the symmetry g± <-> —g±. These two su(2) invariant theories are quite different. The theory 
with <7|| = — g± < is a massive theory with UV fixed point at g\\ = g± = 0, and corresponds to the su(2) invariant 
limit of the usual massive sine-Gordon theory with an empty IR fixed point. The theory with g^ = g± > on the 
other hand is a massless theory with a non-trivial IR fixed point at <7y = g± = 0. In principle the UV behavior of 
the latter theory is not dictated by the action eq. JTSJ; the action is rather an effective low-energy theory. However 
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it is well understood that this flow can arise from the 0(3) sigma model at 8 = ir, which is defined on all scales, in 
particular the UV. A masslcss S-matrix for this flow was given in ref. 20] . These observations lead us to consider 
two different S-matrix theories for the region h ^ 0. The first is the massive theory proposed in ref. Q which when 
h = is the same as the sine-Gordon theory on the line gu = —g± < 0. The other theory, not considered in is a 
massless theory which corresponds to the 0(3) sigma model with 9 = 7r in the limit h = 0. As argued in section II, 
the massless theory is perhaps more consistent with the RG flow since it has the signature of a RG limit cycle on all 
length scales. 



A. Higher k 

In ref. the beta-functions are given for arbitrary level k. The fc-dependence is quite simple and can be scaled 
out of the RG flows by defining a new couplings g = kg and a rescaled RG time I — l/k. In terms of g and I the flow 
equations are the same as for k = 1. Thus the period in I is 2irk/h, where h is defined by Q = — h 2 /16, with Q the 
same as in eq. I|19|) with g replaced by g. However since the one- loop beta function does not depend on fc, it is more 
sensible to normalize the RG invariant as : 

9u - 9 2 \ h 2 

- - (22) 



(l + fc 5 ||) 2 (l-fc 2 gi) 16 



since then Q = g 2 — g\ + 0(g 3 ) and the leading one-loop contribution to Q is the same for all k. One sees that 

h = kh, and thus the period in I is the same as for k = 1: A = 2ir/h when h is defined in the manner eq. Ij22(l . 
In the next section we will propose an S-matrix for the higher k theories. 



IV. THE S MATRICES 



A. Massive case 

In reference Q it was conjectured that the cyclic regime of the KT flows described in the last section has a spectrum 
consisting of a massive soliton and anti-soliton with topological charge ±1. A factorisable S-matrix was proposed 
based on the U q (sl(2)) quantum affine symmetry of the underlying Hamiltonian where quantum parameter q is related 
to the period 2n/h of the RG cycles by the equation)^ 

q = -e-^ 2 (23) 

This value of q is rather unusual for a relativistic factorized S-matrix theory. The other previously known models 
with S-matrices related to a quantum affine symmetry all correspond to q a pure phase, and are models with fixed 
points corresponding to free bosons or their quantum group reductions to minimal models when q is a root of unity. 
S-matrices with a real q however do appear in the study of the XXZ spin chain in the massive antiferromagnetic 
regime, where the solitons are understood as massive spinons. The relation of our model to the XXZ spin chain was 
discussed in Q. There it was shown that in the low energy limit of small momentum, and as h tends to zero, our 
S-matrix agrees with that of the spinons of the XXZ spin chain. However, in contrast to the spin chain, our dispersion 
relation is relativistic on all energy scales, as appropriate to the relativistic field theory of the last section, and this is 
the main difference between the models. As we will see, the interesting properties of the finite-size effects which are 
a consequence of the cyclic RG only appear in the ultraviolet limit for this massive case, and are thus unlikely to be 
properties of the spin-chain since the latter has an explicit lattice cut-off which leads to a non-relativistic dispersion 
relation. 

We now describe in detail the S-matrix proposed in ref. Q. As usual, we let fi parameterize the energy and 
momentum of the relativistic particles as in eq. Q), where m is the mass of the soliton and anti-soliton. To describe 
their scattering amplitudes let us introduce the creation operators A a ((3) for the solitons/antisolitons, where a = ± 
denotes the topological charge. The S-matrix may be viewed as encoding the exchange relation of these operators: 



(24) 
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The S-matrix can be presented in the following form: 



where q is given in eq. (|23[) . and 
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o 





(25) 



C = e" 



-i[5h/2 



({3 = /3i — /?2-) The overall scalar factor is 



2ig 



~ (l-g 4 + 4 "C' 2 )(l~g 2+4n C 2 ) 

1 - g 2 C 2 (l-g 4 + 4 "C 2 )(l-g 2+4 "C^ 2 ) 



(26) 



(27) 



This infinite product is convergent since |g| < 1. 

The above S-matrix satisfies the necessary constraints of crossing-symmetry, unitarity, real-analyticity and the 
Yang-Baxter equation. The factor p can be understood as the minimal solution to crossing and unitarity when q is 
real. 

The arguments leading the above S-matrix are entirely independent of the RG analysis. Indeed the S-matrix depends 
on the couplings g±, gu only through the RG invariant h. The physical mass m also depends on the couplings, though 
in an unknown, non-perturbative fashion. It is therefore a non-trivial consistency check and confirmation of the cyclic 
RG if the S-matrix possesses the cyclicities discussed in section II. Indeed, it does possess the cyclicity eq. © with 
exactly the period 2A predicted by the RG analysis based on the all-orders beta-function. This is easily verified since 
C(/3 + 2A) = C(/3). Recall that the period of the RG at one-loop is h/n which is A/2 3]. Thus, that the all-orders 
beta-function correctly predicts the RG period is an indirect check of it. 

S-matrices with real periodicities in the rapidity have previously been studied in |22l |23|. Indeed, it was realized 
in p2| that such a periodicity may imply an RG limit cycle; however at the time it was unknown what field theory 
the S-matrix described so it wasn't possible to do an RG analysis. In these works, the S-matrix is expressed in terms 
of elliptic functions, whereas our S-matrix is trigonometric. For a single particle with diagonal scattering, crossing 
symmetry and unitarity lead to an imaginary 2iri periodicity, so that an additional real periodicity in j3 naturally 
leads to the doubly periodic elliptic functions. Our S-matrix on the other hand is non-diagonal, and thus is not 2iri 
periodic, which explains why it can be expressed in terms of trigonometric functions. 

One problem with the above massive S-matrix has already been alluded to. As discussed in section II, the above 
periodicity in rapidity implies an RG limit cycle only in the deep UV, which contradicts the flows described in section 
III which have limit cycles on all length scales. As we will see, other problems are encountered when we investigate 
c G ff- We will argue that the massless S-matrix theory on the other hand doesn't suffer from these problems. 



For later purposes it is convenient to express the component Sff as follows, 



ilogS++(/3)=7T+^ 



0h 

Y 



2 sm(n(3h) 



^ n 1 + e nvh 

n— 1 



In the limit where q — > — 1, i.e. /i — s- 0, it is easy to show that l|28|l becomes the 
at the SU(2) invariant point (SU(2) Thirring model): 



(28) 



entry of the sine-Gordon model 



nog s+i (/?) 



7r + f — sin(flx) 

Jo X 



cosh(7rx/2) 



(29) 



B. Massless case 

An S-matrix description of theories with IR fixed points was given by Zamolodchikov and Zamolodchikov |2fl \2l\ . 
An essential ingredient of their formulation are formal S-matrices for only left-movers or only right-movers, denoted 
Sll and Srr. These S-matrices are formal S-matrices for a scale-invariant theory, and encode the information, such 
as the Virasoro central charge, of the infra-red fixed point. When left-right scattering Srl is non-trivial, the UV 
properties are affected. For the known models this collection of S-matrices corresponds to theories with non-trivial 
fixed points in both the UV and IR. 
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We now consider a massless S-matrix description of the limit cycle theory described in section III. Let us param- 
eterize the massless energy momentum for left and right movers as in eq. JJJ. Requiring the two-particle S-matrices 
to correspond to the 0(3) sigma model at = ir in the su(2) invariant limit h —* leads to the obvious proposal: 



Srr(P) 


= S h ((3), 


P 


= Pri 


-PR2 


Sll(P) 


= S h (P), 


p 


= Pli 


-PL2 


Srl(P) 


= S h (l3), 


p 


= Pm 


-PL2 



where all the S-matrices Sh on the right hand side are the same as in eq. I|25l) . 

The above S-matrix again has the periodicity anticipated in eq. In this massless situation, this is a signature 
of cyclic RG on all length scales, entirely consistent with the RG flows. 



C. Higher k 

In this subsection we propose S-matrices for the higher level k models. In the limit where h = 0, the model is 
the well-known isotropic current-current perturbation. In the massive case (gu = — g± < 0) the particles are known 
to have RSOS kink quantum numbers in addition to the spin 1/2 su(2) quantum numbers. The S-matrix factorizes 

= S su{2) <g> Sj$, 0s where S su{ 2) is the su{2) invariant S-matrix of the su(2) Thirring model and S^ os is the 
k-th RSOS S-matrix which is a quantum group restriction of the usual sine-Gordon S-matrix (restricted sine-Gordon 
model) and occurs in the perturbations of the k-th c < 1 minimal conformal series pi) . This factorization is a 
consequence of two commuting symmetries, the su{2) and the (fractional) supersymmetry |2l| . 

Based on the above h = result, it is straightforward to propose an S-matrix for the k-th massive cyclic the- 
ory. When h ^ 0, the su(2) symmetry becomes the U q (sl(2)\ symmetry with the same as above, the fractional 
supersymmetry is unaffected, and both symmetries continue to commute. The S-matrix is then: 

= S h ® S^ sos (31) 

where Sh is the S-matrix described in eq. H25|) . and S^g OS is independent of h. 



V. FINITE-SIZE EFFECTS: MASSIVE CASE 



A. Ground state scaling function 

As discussed in section II, let us place our model on a cylinder of finite circumference R and length L going to 
infinity, and consider the ground state energy E(R). In this section we consider only the massive case. 

Given the S-matrix one can derive integral equations for c c g by means of the thermodynamic Bethe ansatz (TBA) 
[l3| . The derivation of the TBA equations for diagonal S-matrices given in makes no reference to the usual Bethe 
ansatz for the eigenstates of a lattice regularization of the field theory, and relies only on knowledge of the S-matrix. 
Unfortunately for non-diagonal S-matrices the TBA normally leads to an infinite number of coupled integral equations 
for the pseudo-energies. However a dramatic simplification to a single integral equation has been found by Pearce 
and Klumper |24j, and Destri and de Vega [2{|. In particular, Destri and de Vega (DdV) derived equations for the 
usual sine-Gordon theory based on a special scaling limit of the XXZ spin chain. Since as described above it is in 
principle possible to derive the TBA equations directly from the S-matrix alone, wc can take the point of view that 
DdV indirectly have dealt with the mathematical complications of the non-diagonal nature of the scattering in the 
sine-Gordon theory. Since our S-matrix is essentially a sine-Gordon S-matrix with a real q rather than a phase, we 
can thus obtain the TBA equations of our model from that of the sine-Gordon model equations by simply using our 
expression for the S-matrix. 



B. The DdV equations 



In the approach developed by DdV, the effective central charge has the following description in the usual sine-Gordon 
theory in the regime with only solitons in the spectrum: 

CcS (R) = 3m / dp sinh(/3 + ia) log(l + e lZ{(i+ta) ) (32) 

" J — OO 
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where m is the soliton mass and Z is a solution of the non-linear integral equation: 

y* oo 

Z(/3) = mR sinh /3 + 2 3m / d/3' G(/3 - /?' - ia) log(l + (33) 



The kernel G(f3) depends only on the soliton to soliton S-matrix element: 

G03) = ^^logS++(/3) (34) 

Above, < a < n/2, is a regulator needed for convergence. 

For the usual sine-Gordon model with given in ref. |26| . c e g(mi?) is an uneventful function that smoothly 
interpolates between Cfjy = 1 in the ultra-violet (R — > 0) and Cm = in the infra-red (R — > oo). This reflects the 
existence of a c = 1 UV fixed point, a free massless boson, and an empty IR fixed point. The fact that cm = 
is simply due to the fact that the theory is massive and in the deep IR these massive states decouple. This is in 
accordance with the rough picture that c e s is in a sense a measure of the massless degrees of freedom: in the deep 
UV the massive particles are effectively massless and cuv = 1 , whereas in the deep IR the massive particles decouple 
and cm = 0. We have verified the above statements by solving the integral equation (|33|l numerically for the usual 
sine-Gordon model. The UV and IR limits of c c s can also be computed analytically in terms of dilogarithms [TH.!^ . 

Because of the cyclic nature of the RG in our model, we expect a much more interesting behavior in the UV due 
to the non-existence of a fixed point. In the IR on the other hand, the usual argument given above applies since the 
theory is massive and one concludes cm = 0. 

For our model, the kernel has the simple expression: 



1 



4 £t 



cos(n/3/i) 



n—1 



(35) 



which follows from eqs. I)34I28[I . The kernel has the following periodicity: 



G (f 3 + t) = G(/3) (36) 

which is a clear indication of the cyclic RG, with the period exactly as predicted by the RG analysis reviewed in 
section III. It is this cyclicity of the kernel that is ultimately responsible for the non-trivial behavior of c n n in the UV. 

Let us first rewrite the DdV equations so that they resemble the more conventional expressions |13| . and also 
make them more suitable for numerical solution. Since Z{(3) is a real function when f3 is real, let us define the 
pseudo-energies 

e(j3) = -iZ(P + ia), s*{/3) = iZ(fi - ia), (37) 
Then eq. I|33|l can be written compactly as 

e^eo-GQ^L + G^L*, (38) 



where * denotes a convolution 



and we have defined: 



(/*<?)(/?)= / d0 /(/?-/?') <?(/?')• (39) 



The expression for c g becomes 



e (/3) = -imRsmh((3 + ia), (40) 
L[fi) = log(l + e- e <«), (41) 
G (/3) = G(fi), = G(f3 + 2ia), (42) 



6 f°° 

c cB (mR) = — SRe / dp e (P)L(P) (43) 
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We have solved the above equations numerically and verified that c e g is independent of a. A detailed presentation 
of our numerical results will be given below. For the further development of analytic results, the value a = tt/2 is 
preferred since then Eq — mR cosh (3 is the standard pseudoenergy for a free theory, and when the kernel G is zero, 
which occurs at the free fermion point of the usual sine-Gordon model, one has e — Eq and the expression for c e ff is 
the standard free fermion one: 

(\m f? /*°° 

c ite JmR) = —P d[3cosh(3 \og(l + e~ mRcosh(3 ) (44) 

7T 2 J-oo 

In fact, in the usual sine-Gordon theory, c c g is rather well approximated by Cf rcc at all values of the coupling. 
Henceforth a will be set to tt/2 in all analytic computations, and unless otherwise stated: 

e ((3) =mR cosh (3. (45) 

It is convenient to define the function 

a(0) = e(fi) - e (/3). (46) 

Since the kernels Go,i are periodic functions of (3 with period 2ir/h, then so is u(/3), and it can thus be expanded as 
a Fourier series: 

oc 

a(/3) = °m mhfi (47) 



n— — oo 



We can then proceed to rewrite the integral equations as non-linear algebraic equations for the various Fourier 
coefficients. First we expand the kernels: 

h °° 1 

= 2; E 90 > neinh *> 50 - = TTd^ (48) 
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Then the integral equation l|38(l can be written as 

o"„ = -go, n L n + gi t nL- n , (49) 

where we have defined 



L n = — J df3e- inh ?L(l3), 



(50) 



with L(j3) as before but now viewed as a function of a. In deriving eq. 149(1 we have used that (L n )* = L n , which is 
a consequence of (L((3))* = L{—(3). From equation (|4^|) one can derive the following relation: 

9l,n -nirh„ /mi 

er„ = cr_„ = -e cr_„. (51) 

90,n 

which implies the absence of the zero mode, i.e. a = 0. 



C. Exact Numerical results 



The numerical work is simplified by splitting the various functions in the model, i.e. / = e,a,L, and the DdV 
equations into their real and imaginary parts (e.g. f = f + if")- Moreover the symmetry properties satisfied by 
these quantities further reduce the computational resources of the program. The kernel G(/3) in eq. (|35[l is an even 
function of /3, and Z((3) in eq. (|33f) is a real quantity. These facts imply that Z((3) is an odd function of /3, which 
in turn determines the parity properties of the remaining functions in the theory. Namely, the real parts, /'(/3), are 
even, while the imaginary parts, /"(/3), are odd functions of /3, and therefore /(— j3) = /*(/?). 

The solution of the DdV equations has been found using two different methods: 1) recursive solution of the eq. i|38|) 
for the functions cr'(j3) and u"(/3); and 2) finding the roots of the non-linear system l|49|l for the modes <r n subject to 
the constraint l|51|l . The results obtained by these two methods are in full agreement. 
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In fig^we show the solution to eq. I|38|) in the case where h = 1 and several values of mR ranging from mR = 10~ 2 
to 10 -8 . These results confirm that a' is an even and periodic function of with period 2ir/h, while a" is an odd 
periodic function with the same period. We also observe that the amplitude of the oscillations is rather small, of order 
10~ 2 in this case, which is a general feature except for special regions of mR where the amplitudes are of order 1. 

Fig-Elshows the values of c e g(mR) corresponding to h = 0.5, 1, 2 in the range mR = 1 — 10~ 14 . In the IR region 
c e ff goes to exponentially fast and independently of the value of h, as expected for a massive model. As soon as 
mR < 10 -2 , c c ff starts to oscillate around c = 1 with a period n/h in the scale, i.e., 

c(mR) ~ c(e* lh mR) (52) 

Indeed, as seen in the figure, a complete cycle for the h = 0.5 model contains two and four cycles of the h = 1 and 
h = 2 models respectively. Another feature of c e g is that the amplitudes of the oscillations increase to-wards the UV, 
which is especially pronounced for h = 0.5, and this is why eq. I|52|l is not an equality. The RG analysis performed 
in section III suggests a periodic c c g (mR) function with period 2ir/h in the scale, and this is twice what is observed 
numerically. This will be explained analytically in the next subsection. 



D. Approximate analytic solution 

In this subsection we derive an approximate analytic expression for c e ff, explaining the properties observed nu- 
merically. Later we will show that it agrees remarkably well with the exact numerical solution to the full non-linear 
problem, and arguments for why this is so will be given. 

We first make a linear approximation to the integral equation (14911 . Using 

L(f3) = Lo(f3)-K(P)a(f3)+0(<T 2 ) (53) 

where 

L ((3) = log(l + e-^)), K{J3) = 1 + ( l o(/3) (54) 

then to order a: 

L n = £o,n ~ ^ K n - m a m (55) 

m 

where Lo,n and K n are Fourier transforms of Lq and K: 

L 0>n = — / df3 e-^LoiP), K n = — dp e~ inh ^K{l3) (56) 

Using eq. (|55|l in eq. 149(1 one finds the linear equations: 

^ A n . m <j m = B n (57) 



where 



In this approximation, 



An : m — &n,m Q0,7i-^n—m "T" QX^n-^—n—m (58) 
B n = gi tn Lo,-n — 30,raio,n 



22 

Ceff (mi?) sa c llcc (mR) } o n K' n (59) 



where K' n is the following Fourier transform: 



K = ^ / d(3 e-WK'ifi), K'(fi) = e ([3)K(P) (60) 
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All of the needed Fourier transforms Lo,n, K n , K' n are computed in appendix A. There the mR dependence is 
expressed in terms of the logarithmic variable 

s = log (2n/mR) (61) 

In the deep UV where R goes to 0, s is very large. 

From the results in appendix A one finds in the UV limit of large s: 

K « A( 8 _ 7 _iog 2 ) (62) 

Z7T 



where 7 is Euler's constant and 



K ^ w . /0 . jfte(e*" fts Ci(l-^)) 

7rcosh(n7r/i/2) v ' 

d(z)^(l-2- z )C(z) (63) 

with £ Riemann's zeta-function. Thus Kq is large compared to K n ^Q and this justifies the further approximation of 
keeping only the K terms in eq. <|571) . (if„^o is further suppressed by e ~ n7Th / 2 ^ Using eq. (|5T)l . this essentially 
diagonalizes eq. (f5*7|) and one finds 



B„ 



(64) 



1 - 2go,nK 

Using eq. (|64(l and the results in appendix A in eq. (|59J) one finds 

CcS (mR) = c lrcc (mR) + - £ J fe ^ ^— + 0(e- 2s ) (65) 

which is meant to be valid in the UV. The contribution Cf ree rapidly approaches 1 in the UV. 

A direct consequence of this equation, together with 1)61(1 . is that the amplitudes of the oscillations of c e s(mR) 
increases to-wards the UV because of the presence of a pole at the critical distance i? Ci „ 

mR c , n = Tre" 7 exp (- ~ (1 + e™ 7 ^)) (66) 

One sees that c s{mR) becomes singular at i? Cj „ in this linear approximation. The double exponential in (|66|l implies 
that these scales are extremely small. For example the highest value of R Ct i(h) is attained at h = 0.407 where 
mi? Cj i = 7.14 10~ 16 . Therefore at this scale we expect our approximation to break down. In the next subsection we 
shall study in more detail these critical points. 

In fig. |3 we compare the result (|65|l with the exact numerical one for h = 0.8 in the range mR = 1 — 10~ 20 , 
observing a very good agreement over a wide range of distances. In this case the first critical distance is given by 
mi? Ci i = 3.06 10 -23 . Eq. (|65)1 explains the oscillatory behavior of c e s(mR) with a period ir/h in the scaling variable s. 
The lack of periodicity is due to the denominator of this expression which becomes more important as one approaches 
the critical points R c , n - 

The reason why the linear approximation is so good is that the higher order terms in the expansion of L are in fact 
strongly suppressed in the deep UV (see Appendix B for details). 

In fig.0|we give a similar comparison for h = 0.4 over a range that contains two critical distances mi? Cj i = 7.09 10 -16 
and mR c ,2 = 5.33 10" 46 . 

Finally consider the h — * limit, where the period A becomes infinitely large. The singularities at i? Ci „ are pushed 
deeper and deeper to the UV and approach zero. One thus approaches a smooth flow between cjjv — 1 and cir , = 
as expected for the isotropic limit of the massive sine- Gordon theory. 

A curious property of ea. H65l) is that for the special value h — hi = 2ir/ log 2 the linear contribution to c c ff vanishes 
as a consequence of the identity ^i(inhi) = 0, Vn. Thus the value of c e g is given simply by Cf rcc . The exact numerical 
solution also satisfies this result to order 10 -16 which was our machine precision. For h nearby hi the value of c c g — Cf ree 
is small but never zero. This fact is a consequence of the "absence" of the prime number 2 in Euler's product formula 
for the zeta function Ci( z )- For values of h 7^ hi the corrections to Cf rcc in eq. (|65|l are never zero since C,(ih) does 
not vanish. The later quantity is related by the Riemann's functional equation to £(1 — ih) which does not vanish as 
well. The latter result was used by Hadamard and de la Vallee-Poussin to prove the prime number theorem |28j . 

In summary, in the massive case the effective central charge c c ff exhibits nearly periodic structures over a large 
range of scales in the ultra-violet, with a period consistent with the RG predictions. However as one approaches the 
scales R Cin in the very deep ultraviolet, c g develops a singularity and the periodicity eq. I|14|l is violated. In the 
analytic approximation this singularity appears as a pole in s in the equation 1)65(1. These features suggest that the 
massive theory is not entirely consistent with the RG flows, as alluded to in section IV. 
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Exact result 
— Linear approximation 
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FIG. 3: Comparison of c c s(mR) computed with the analytic approximation eg, 165H and the exact numerical results for h = 0.8 




Exact result 
— Linear approximation 
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FIG. 4: Comparison of c c a(mR) computed with the analytic approximation ea. l|65|l and the exact numerical results for h = 0.4 
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FIG. 5: Conformal map L(e) = log(l + e e ) from the cylinder with a cut at e' < and e" = (2p + l)n into the strip 
— 7r < 3mL(e) < n. 



E. Coexistence regions 

The approximate analytic solution of the DdV equations presented above explains the main features of the behavior 
of c e ff, namely its periodicity and scale dependent amplitude. For most values of mR the mode amplitudes a n are 
rather small for all n, which explains the high accuracy of the analytic approximation. However at the critical values 
R c .n the n th mode blows up and the analytic approximation is no longer valid, meanwhile the numerical solution still 
gives finite results. 

Before we focus on this problem let us consider an important technical point regarding the definition of the logarithm 
L(/3) in (|4*T1) . The imaginary part of L{[3) is the phase of w(/3) = 1 + er e ^\ For — oo < /3 < oo the variable w(f3) 
describes a trajectory in the complex plane, which may pass through the cut of the logarithm, changing branches. 
As long as the modes a n remain small enough this never happens. We have found that the solutions of the DdV 
eqs. belong to the principal branch of the logarithm, i.e. — n < ^smL(/3) < n. Under this condition L(e) becomes a 
conformal map from the cylinder with a cut at e' < and e" = (2p+ into the strip — 7r < 3raL(£) < w, as shown 
in fig. [5] The cut corresponding to p = translates into, 

oo 

e' = £q + a' n cos(nhf3) < 0, a' n = — t&nh(nha)cr", 

n=l 

oo 

e" = 4' + ^<'sin(n/^) = tt, (67) 

n=l 

and for a = ?t/2 (£q = 0) these equations simplify. Since just one mode dominates the behavior of the system near 
the corresponding singularity, let us consider the single mode equations, 

e' = mR cosh (3 + a' n cos(nh(3) < 0, a' n = — tanh(n/i7r/2)cr", 

e" = a'^sin(nh/3) = tt, (68) 

which implies that < ir. Now we can solve the eqs. (|49|l for a single mode and search for the values of R around 
R c .n satisfying cr" = ±tt. The results are shown in fig. 

This study of the logarithm solves the problem of the unbounded growth of a'^ and consequently c c e(mR). In 
the IR there is a unique solution which evolves continuously to-wards the UV until it hits the cut, i.e. |cr"| = tt, 
beyond which there is no solution. The boundary of the "IR phase" has a sawtooth shape which on average follows 
the analytic estimate (|66|l (see Fig. |BJ. For smaller values of R < i? Cj i one finds again a unique UV-solution which 
can be smoothly continued to-wards the IR until it hits the cut again. The boundary of this UV-solution has also a 
characteristic sawtooth shape as shown in Fig. |SJ The IR and UV phases may coexist, meaning that there are two 
solutions for the same values of mR and h. The regions with no IR/UV coexistence may contain isolated solutions 
non connected to the IR or UV regions, or they may not contain solutions at all. We have encountered both kinds 
of situations. The previous discussion concerns the truncation of the DdV equation to a single mode. Including all 
the modes yields an infinite number of coexistence regions, one for every mode, as shown already by the analytic 
approximation. The whole picture is as follows. Moving from the IR to-wards the UV, the first mode a'{ passes 
through its critical region and gets damped beyond it. Then, the second mode takes over and after a crossover region 
the period of c c s doubles. This second mode also reaches its critical region and the process repeats itself with the 
third mode ( see fig. 0. Once more the analytic result fits amazingly well the exact numerical result including the 
proximity of the coexistence regions (for clarity we have truncated the analytic part very close to the pole). 
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FIG. 6: Coexistence of the IR and UV phases. 




mR mR 



FIG. 7: c a g(mR) and the a[ and a' 2 in the region of the "singularity" for h = 0.4. 
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Finally, fig. [7| presents a comparison of the exact solution and the single mode approximation near the first singular 
region. At h — 0.4 there is IR/UV coexistence as can be checked in fig. [5] Hence the IR and UV curves c c ff(mR) 
overlap around mR ~ 10~ 15 (fig. EH). Fig. [3d shows that the first mode of the exact solution, er", dominates over the 
second mode, a'^, which in the UV becomes more relevant. Notice also that a" is very well described by the single 
mode approximation. 



VI. REMARKS ON c off FOR THE MASSLESS CASE 



The structure of the TBA equations for massless theories has been understood in a number of models [T^ . I20I [27J • 
Unfortunately the massless DdV form of the equations are not known for the 0(3) sigma model at 6 = tt, thus unlike 
the massive case, it is unknown what the h = limit of the TBA equations based on the S-matrices in eq. H30|) should 
be. In spite of this, based on the known structure of massless TBA equations we can anticipate the general properties 
c e ff for theories based on massless S-matrices with limit cycles. 

When the LR scattering is trivial, the S-matrices Sll and Srr describe a conformally invariant theory. Let us 
introduce pseudo-energies Sl, En separately for the left and right movers. As for other models, when the LR scattering 
is trivial the TBA equations are two decoupled integral equations for e L and £ R , each equation having the same form 
as a massive theory, the only difference being in the source terms £q,£ R which now reflect the massless dispersion 
relations. Thus, when the left-right scattering is turned off, the TBA equations should just be two decoupled DdV 
equations of the kind described in section V. Turning on the left-right scattering Slr then couples the L,R integral 
equations. 

Consider a toy model with only one particle with Sll = Srr = 1 but a non-trivial Slr, and define as usual the 
kernel as 

G((3) = ^dplogS LR (69) 

Z7TI 

An example is the flow between the tricritical and critical Ising model p^ |. Define the free pseudoenergies: 



mR 



em 



mR 



The structure of the TBA is the following: 



where 



The effective central charge is: 



e L = Eq — G * L R 
s R = s R -G*L L 



L,R 



log (1 + exp(-e L ' R )) 



C e ff (mi?) 



d/3 



(e R L R 



e L n L L \ 



(70) 

(71) 
(72) 

(73) 



When the left-right pseudo-energies are decoupled, then by simple shifts of the rapidity one can show that c e g(e l R) = 
c c ff(R) for any / indicating a scale invariant conformal field theory with c e g independent of R. When there is non- 
trivial LR scattering, which leads to the terms in eq. (|71|l which couple L and R, then c c ff becomes a function of 
mR. 

Now consider an S-matrix that is periodic in rapidity, as for our model, such that G{(3 -I- A) = G(f3), where as usual 
A is the period of the RG. An example of such an S-matrix for a single scalar particle was given in [23. The due to 
the periodicity of the kernels G, c e g is now invariant under a discrete scale transformation e x . To see this, one notes, 
using the periodicity of the kernels, that Sl((3— A) and er((3 + A) satisfy the same integral equations as £l,r but with 
R replaced by e x R. Namely, 



e L ((3 - A, R) = e L ((3, e x R), e R ((3 + A, R) = s R (j3, e x R) 
Using these equations in eq. i|73f) one can verify the periodicity of c e s: 

c cS (e x R) = c Gff ( J R) 



(74) 



(75) 



The above exact periodicity of c e ff over all length scales is rather remarkable and entirely consistent with the RG 
argument that led to eq. (|14|) . Also, it indicates that unlike the massive case, there is no scale dependence of the 
amplitudes of oscillation of c c s and thus no singularities of the kind encountered in the last section. 
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VII. CONCLUSIONS 



We have analyzed the finite size effects, equivalently the finite temperature effects, in the ground state energy or 
effective central charge c c ff for the massive S-matrix theory proposed in The massive S-matrix was conjectured 
to describe the cyclic regime of the Kosterlitz-Thouless flows The S-matrices are periodic in the rapidities, with 
a period consistent with the RG analysis. This periodicity leads to novel finite size effects which also have periodic 
structures consistent with the cyclic RG. We obtained an approximate analytic expression for c c g expressed in terms 
of Ricmann's zeta function eq. (j65JI . which has log-periodic behavior and agrees remarkably well with the exact 
numerical solution. 

For the massive case, there is a wide range of scales in the UV where c e g oscillates. However in the very deep UV 
a singularity is encountered, and in some regions c e s is not well defined (coexistence regions described in section V). 
This suggests the massive S-matrix theory needs an explicit UV cutoff in the very deep UV. For the massive case c e ff 
also lacks a periodicity on all scales as would be predicted by the RG. 

For the massless case the TBA equations are unknown and we could not carry out the comparable analysis. However 
we argued on general grounds that the massless case should have some appealing properties. The c e g should be exactly 
periodic on all scales, consistent with the RG. We hope to report on this in a future publication. 

The log-periodic dependence of c c g(mR) on mR is reminiscent of the log-periodic fractal dimensions of self-similar 
fractals, such as the Cantor's triadic set. Self-similarity amounts to discrete scale invariance, which in turn leads to 
log-periodic behavior in a wide variety of problems such as DLA (diffusion-limited-aggregation clusters), earthquakes, 
financial crashes, etc. [2(1 IsnHsij l . 

On a broader note, there appears to be a network of deeply interrelated concepts and techniques, namely RG 
limit cycles, discrete scale invariance, complex exponents, fractals, log-periodicity, quantum groups (with real q), zeta 
function regularizations, number theory, etc, whose full significance needs to be clarified. 
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Appendix A 

In this appendix we compute the coefficients F n = {Lo. n , K n , K' n } of the functions F((3) — {Lq(P), K((3), K'(f3)} 
needed in section V-D. They are given by the formula 

Fn = ^~F{nh) (76) 

where F(t) is the Fourier transform of F: 

/>oo 

F(t) = / d/3 F{(3) e- ltlB (77) 



A. Fourier transform of K (J3) 

The function K((3) is given by 

K{P) = l + exp(rcosh(/?)) (78) 
where r = mR is the scaling variable. Since K{(3) is an even function we can write its Fourier transform l|77|) as 

^ f°° p^P 

K{t)= / d/3- — — (79) 

oo l + exp(rcosh(/3)) 
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For t > one can evaluate (|79|l as the sum of the residues of the poles of K(f3) in the upper half plane 3m/3 > 0. 
These poles are the zeros of the denominator of K((3) and are given by 



d%, = ±e p + i-(i + 2 P '), P ,p'>o 



(80) 



sinh 



The residue formula yields 



K(t) 



1 /7r(l + 2p) 



E 



2n 



cos(t8 p ) 



rcosh(7rf/2) ^—^ cosh( 



(81) 



For t > and real this scries diverges. We shall regularize it using the Riemann zeta function C,{z). For simplicity let 
us consider the UV limit r — > where P has the expansion, 



9 P = s + log(l + 2p) 



+ 0{e- 4s ) 



where s — log(27r/r). The series (|81|l becomes 

1 



K(t) 



E 



COSh(7rf/2) ^ 



(1 + 2p) 2 

„- 2 . ( 2 - «) 



,(l + 2p) 1 - it (l + 2p) 3 -" 

The sums in p can be done using the Riemann zeta function £(z) 

«*) - E i 



+ (t - -t) 



+ 0(e" 4s ) 



71=1 



together with the relation: 



oo 1 

y L 

^ (1 + 2 



(1-2- 2 ) CW^CiW 



This gives 



K(t) = 



1 



cosh(?ri/2) 



- (1 + 2 P y 

[e lts [Ci(l - it) - e- 2s {2 - zt)Ci(3 - it)] + (t -> -t)] + 0{e' 4s ) 



(82) 



(83) 



(84) 



(85) 



(86) 



Strictly speaking eas. (|84|l and (|85|1 are only defined when die z > 1. However the Riemann zeta function is defined 
in the entire complex plane except at z — 1 where it has a pole. Our use of £(z) is completely analogous to other 
situations in quantum field theory where it serves to regularize formal expressions. 

Equation (|75|l gives the Fourier coefficients 



2tt 



K(nh) 



(87) 



In the UV limit the terms proportional to e 2ps (p = 1,2,.. .) vanish very rapidly and one is left with the oscillating 
terms that depend on d(l ± it). The value K Q is special since £(z) behaves near z = 1 as 



C(z) 



z- 1 



+ 7 



where 7 is the Euler's constant. Using this equation the value of Kq is given by 



K = — ( s-7-log2 



7C(3) 



-2« 



+ 0(e" 4s ) 



(88) 



(89) 



The linear dependence of Kq with s in the UV limit can be easily understood using the "kink method" proposed 



suy v. 

by Zamolodchikov to study the UV limit of usual massive models [13. From ea. H77l76[> we have 



K "=2^ 



dd 



1 



1 + exp(r cosh/3) 



(90) 
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which in the limit r <C 1 can be approximated as 

h f°° 1 

A 'o = - / d $ TT n^7\ ( 91 ) 

7T J 1 + exp(eP s ) 

where s' = log(2/r) 3> 1. From this expression it is clear that the integrand remains close to 1/2 for < (3 < s' and 
that it vanishes exponentially for /3 3> s'. Hence ~ hs'/(2w) which agrees with ea. 1)89(1 up to constant terms since 

s' = S — log7T. 

B. Fourier transform of K' 

The function K'(j3) is defined as 

K'((3)=e (m(P) (92) 
Following the same steps as in the previous subsection we find the Fourier transform, 

2n °° 

*' {t) = S inh(7rt/2) ^ tanh(0p) S[n(t0p) (93) 
In the UV limit this can be approximated as 

K'(t) = - g . nh ^ /2) [e lts [Ci(-it) - e" 2s (2 - »t)d(2 - it)] - (t -> -t)] + 0( e - 4s ) (94) 

where we have included explicitly the term of order e _2s as in ca. H86|) . The value of this function at t = is given by 

F'(Q) = -4C(0) log 2 = 2 log 2 (95) 

where we have used that £(0) = —1/2. 

Using the kink method one can see that this result is equivalent to the following integral 

1 f°° e 13 

21og2 = - / d/3 - 1 r T - (96) 

6 2 j— oo 1 + exp(e^/4) 1 ' 

which can be checked numerically to any precision. 

C. Fourier transform of Lq{0) 

The function Lq{0) is defined as 

L (/3)=log(l + e -'' cosh ' 3 ) (97) 
To compute its Fourier transform we note that 

-r±L ((3) = ±L (l3)=K'(l3) (98) 
Hence the Fourier transform of Lq{(3) can be obtained integrating ea. (|98|l in s, obtaining: 



L (t) = -- 



sinh(7ri/2) 



e Its ( ~G (-it) + ie" 2s Ci (2 - it) J + e~ iu ( j(i (it) - ie- 2 Xi (2 + it) 



(99) 



t v J \t 

The integration constant is zero as we shall next show. Let us consider the value of Lo(0): 

Lo(Q) = 2slog2-log2 log(27r 2 ) (100) 

where we have used that C'(Q) = ~ \ log(27r). Notice that the leading term in s coincides with the value of -FT'(O) given 
in ea. l|95|) . but there is a constant term that appears in the limit t — > in (|99|l . Using the kink method eq. I|99|) 
amounts to, 



lim 

r— >0 



/>oo 

2 d/91og(l + exp(-rcosh^)) - 2s log2 
Jo 



= - log2 log(27r 2 ), s = log(27r/r) (101) 



which we checked numerically. This shows that the integration constant relating K'(t) and Lo(t) vanishes. 
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Appendix B 

In this appendix we justify the analytic approximation made in section V-D of the non linear ea. (|49|l . 
The full expansion of the function L(f3) is given by 

L(/3) = L (P)-J2 Z T K^(f3) (102) 

a—l 

where 

K W(B) = —, KM =K= (103) 
Using ea. H102|) . the expansion of ea. (|49|l to all orders in a becomes, 

^A n>m (T m = B n + ^- C n,E^ "mr-O-m. ( 104 ) 

m a— 2 mi,-,m fl 

where A n ^ m and B n were defined in ca. (|58|l and 

JfW = _ / d(3 e~ mh0 K {a \f3) = —K^ a \nh) 
2?r J_ 00 2vr 

The analytic approximation of section V-D consists in i) eliminating the terms 0(a n ), with n > 2, in ea. (|104|l yielding 
eq.(J5S}, and ii) considering only the coefficient Kq which depends linearly on the scale s (recall ea. (!62|l % l. The latter 
asumption is justified if the cr„'s are small enough. As for the first assumption we shall show below that the terms in 
eq. (|104f> containing high powers in a n do not depend linearly on the scale s, and can thus be neglected in the UV. 
The Fourier transform of can be found by the techniques of appendix A and it reads, 



£w (t) = _(z!) y e «». r (» + 1 ~ p Cl(0 + 1 _ uu) x ( H c °t ( 1^' ^0,2 

\ -7T / ^ Y{\-itv) v ; \ ^/smh(7rf/2), a = 1,3,... 



i/=±i 



The case a = agrees with eq. (|86fl . Let us focus on the cases where a > 1. Only for t — one may get a linear 
dependence on s, namely 

^)(0)=a!Ci(a+l)x(~^C' fl = • ■ • (1Q7) 

1 4( ¥ ) s, a - 1,3,... 

Returning to ea. H104(l . the only terms that may contribute linearly in s, are 



\- 2h( 1 (a)i a ( ^ a 

/ ^ a+1 I 90,n / J On.rn\-\- ...-\-m a ^mi ■ ■ • @m a / ^ ^-n.mi-j-...+niQ ^mi ■ ■ ■ &n 

> 2, even \ mi,...,m a mi,...,m a 



mi ,. . . ,m a 



(108) 



a>2, even 

Using eq.JSTJ one can show that 

^ &n.m\-\- . . .+m a ®m\ • • • @m a ( 1) ^ ^ & —n.m\-\- . . .+m a &m\ • • • ^m a (^09) 



which together with ea. l|48|) . namely g\ n — e ra <7o,m implies the vanishing of l|108[) . In summary, all the linear 
dependence in s of the eq. (|104fl comes from the coefficient B n . 
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